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Abstract: It is shown that the amplitude for reflection of a Dirac particle with arbitrarily low momentum incident 
on a potential of finite range is —1 and hence the transmission coefficient T = 0 in general. If however the potential 
supports a half-bound state at momentum fc = 0 this result does not hold. In the case of an asymmetric potential the 
transmission coefficient T will be non-zero whilst for a symmetric potential T = 1. Therefore in some circumstances 
a Dirac particle of arbitrarily small momentum can tunnel without reflection through a potential barrier. 


INTRODUCTION 

The results for scattering at arbitrary low energy E in one dimension in the Schrodinger equation are well known. 
If the potential V (x) is sufficiently well behaved at infinity, then the reflection coefficient at zero energy is unity and 
the transmission coefficient is zero [Q unless the potential supports a zero energy resonance (a half-bound state). In 
that case the transmission coefficient is unity and there is no reflection provided that the potential is symmetic. Bohm 
calls this a transmission resonance Q. These results have been generalised to asymmetric potentials |^, [Q. In this 
paper we repeat the analysis for the Dirac equation.To some extent this has already been done by Clemence in 
the mathematical literature in his analysis of the Levinson theorem but we approach the problem as physicists. Our 
results show that transmission resonances will occur in the Dirac equation even for the case where the potential V (x) 
is everywhere positive and thus represents a potential barrier. 

The potentials V(x) we shall consider are smooth and of finite range. In non-relativistic systems for such potentials, 
scattering states with continuum wave functions have E ^ 0 whereas bound states with normalisable wave functions 
have E < 0. A half-bound state or zero energy resonance in non-relativistic scattering occurs when the potential 
supports a bound state of energy E = 12711 in the limit k —> 0: the corresponding wave function thus becomes a 

continuum wave function. An example of this is when a square well is sufficiently deep to just support the hrst odd 
bound state: the resulting wave function describes a non-normalisable half-bound state which corresponds both to a 
particle of arbitrarily low energy incident on the potential from the left and also to a particle of arbitarily low energy 
incident from the right. 

In the relativistic Dirac equation, the notion of half-bound states is more subtle. For a free Dirac particle, there 
exists a gap E < \m\ which separates the positive and negative energy continuum states: the positive energy states 
correspond to particle states and the absence of negative energy states (hole states) describe anti-particles. On the 
introduction of a potential V(x) this gap becomes distorted and bound states now occur between E = —m and 
E = m. A potential which is attractive to particles and supports a half-bound state at E = —m or a potential which 
is attractive to anti-particles and supports a half-bound state at if = m is called a supercritical potential: thus the 
Dirac equation has half-bound states at both E = —m and if = to in contrast to the Schrodinger equation where 
these only exist at if = 0. It follows also that we should talk of zero momentum resonances in the relativistic case 
rather than zero energy resonances. 

In the following sections we discuss the one-dimensional Dirac equation using a two-component approach and 
establish the formalism needed for the consideration of scattering and bound states. We will then prove that Dirac 
particles with energy if > to and arbitrarily small momentum incident on a potential of finite range will be completely 
reflected unless the lower component of a particular wave function vanishes. If the potential supports a half-bound 
state at the threshold energy if = to this condition is shown to be satisfied. In this case there will be a non-zero 
transmission coefficient in general whilst for a symmetric potential, there will be a transmission resonance: the particle 
will tunnel without reflection. In particular, we confirm our previous result 0 that solutions of the Dirac equation 
exist in which a particle of arbitrarily small momentum can tunnel completely through a potential barrier. In the 
Appendix we illustrate our results by considering an asymmetric potential which is soluble analytically. 


THE TWO-COMPONENT APPROACH 

Following an earlier paper Q we take the gamma matrices jx and 70 to be the Pauli matrices ax and az respectively. 
Then the Dirac equation for scattering of a particle of energy E and momentum k by the potential V(x) is 
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to obtain the coupled differential equations 
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For a free Dirac particle of momentum k the solution is 


A 

B 


^ikx ^2 _ ^2 _ ^2 


( 2 ) 


(3a) 

(3b) 


A = 


ik 


E — m 


B = i 


E- 


m 


E — m 


B = 


(^) 


B 


( 4 ) 


Suitable choices for A and B will facilitate future calculations. For threshold problems where E ^ m, choosing 
B = —ik leads to A = E + m and the free particle wave function ip can be written apart from a normalisation factor 
as 


It is clear that in this form the top and bottom components do not simultaneously tend to zero as if —*■ m, fc —> 0. If 
on the other hand we were interested in threshold wave functions where E —>■ —m then choosing B = E — m leads to 
A = ik and the free wave function can now be written (again up to normalisation) as 

’>=(e‘-™)'“' < 6 > 


S-MATRIX FORMALISM FOR THE ONE-DIMENSIONAL DIRAC EQUATION 

The S-matrix formalism for scattering in one dimension for the Schrodinger equation is well known and covered in 
a large number of texts (e.g [|^, ^). The same arguments are applicable for the Dirac equation in one dimension 
1 ^ and here we will summarise a number of the more important results in the context of a relativistic equation. 

We adopt the usual formalism for a Dirac particle incident from the left scattering off the piecewise continuous 
potential Vix) of finite range where U = 0 for |a:| > f where the asymptotic solution ipi{x) of Eqs. (3) for particles 
incident from the left with momentum k and energy E using Eq. is 

which defines the (left) reflection amplitude l{k). We can also define the (left) transmission amplitude ti{k) 

ipi ^ tiik) , x^oo ( 8 ) 

We can similarly define the asymptotic wave function for particles incident from the right as: 

Ipr trik) ^ 


—ikx 


— 00 


(9) 


2 






thus defining the right refiection and transmission coefficients r(k),tr{k). 

The left scattering coefficients and the right coefficients can be simplified further. If we had two independent 
solutions of the Dirac equation 


ipi = 




( 11 ) 


then the Wronskian of the solutions V'l) V'2 of the first order linear differential equations of Eqs. (3) defined as 0: 


W[tpi,tp 2 ]{x) = fi{x)g 2 {x) - f 2 {x)gi{x) 


( 12 ) 


would satisfy W'{x) = 0 with the Wronskian W{x) constant and non-zero. (When fc = 0 it is easy to see that any two 
solutions are not independent and W = 0). We can now evaluate the Wronskian W{^i,'tpr){x) as x ^ ±oo to give 


ti{k) = tr{k) = t{k) 


(13) 


So there is only one transmission coefficient t{k). 

The general solution of the Dirac equation 'ijj{x) can thus be written as a linear combination of ipi and ipr- 

1p = Alpi + B Ipr (14) 

The asymptotic solutions are now found to be 

V- - 41 ™ B ^ ™ ^ (15) 

^ ^ ^ ^ ^ ^ 

where 

A{k)=At{k) + Br{k) , B{k)=Al{k) + Bt{k) (17) 


The coefficients A and B are the amplitudes of the incoming waves for particles arriving from x —oo and x —>■ oo 
respectively. Conversely, the coefficients A and B are the coefficients of the outgoing waves for the transmitted or 
reflected particles. We can now introduce the matrix S{k) which allows us to calculate the outgoing amplitudes in 
terms of the incoming amplitudes. 



S{k) = 


( t{k) r{k) \ 

1 m m) 


(18) 


The flux j Is given by 

j = ilj{x)jxtpix) = iip{x)axtp = i'ip{x)azcrxip = -ip\x)ay'il}{x) 
Using equations (W) and & we consequently find that 


j = 2k{E + m){\A\'^ — \B\'^) x —oo 
j = 2k{E + m)){\A\'^ — IBI"^) x ^ oo 


(19) 


( 20 ) 


The conservation of flux gives us the condition 

|yl|2 + |S|2 = |i|2 + |5|2 


( 21 ) 


Also 

|i|2 + |5|2 ^ (A* B*)S{kAS{k) |AP + \B\^ 
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Hence S{k) is a unitary 2x2 matrix. From equation (p^), this imposes the following conditions on the matrix elements 
of S'(fc): 


T{k) + L{k)=T{k) + R{k) = l (22) 

t{k)r*{k) + t*{k)l{k) = t*{k)r{k) + t{k)r{k) = 0 (23) 

where T{k) = |t(fc)P is the transmission coefficient, L(k) = \l{k)\'^ is the reflection coefficient for a particle incident 

from the left and R{k) = |r(fc)p is the reflection coefficient for a particle incident from the right. It also follows that 

m\ = \rik)\ (24) 

Additionally, if the potential is symmetric, i.e. V{x) = V{—x) (see D below) then 'tp'{x) = azip{—x) is also a solution 
(this is Eq. (^) below). The asymptotic wave function ip'ix) can be found from ^p{x) using equations dl), (0) by 
the substitutions A B and A B. This implies that S{k) must also be symmetric and consequently 

r{k)=l{k). (25) 


in this case. 

The last property we wish to illustrate is the behaviour of the amplitudes t{k), l{k) and r{k) at k = 0. By taking 
the complex conjugate of Eqs. (2-5) with negative momentum —k, we see that ip*^{—k,x) has the same form as 
tpi^rikjX). This in turn implies that 

t*{-k)=t{k) , l*{-k)=l{k) , r*{-k) = r{k) (26) 

So we see from Eq. ( 0 ) that all the amplitudes Z(0), r(0), t(0) are real. This will be of importance for the next 
section. We also have from Eq. ( p^ that 

r(0) = -1(0) or t(0) = 0 (27) 

It follows from Eqs. (0) and (p^ that for symmetric potentials 

r(0) = /(O) =0 or t(0) = 0 (28) 


We discuss this further in D below. 


REFLECTION AND TRANSMISSION PROPERTIES AT ZERO MOMENTUM 

A. The General Case 

Our approach will follow that presented for the Schrodinger equation by Senn ||l^ . When a Dirac particle is incident 
from the left scattering on the potential V{x) of finite range so that V{x) = 0 for |x| > the solution of Eqs. (3) 
in Region I a; < —^ for particles incident from the left with momentum k and energy E is just 


, s , , E + m 

^ =i’i=[ _,k I® 


' ;(fc) ( ^ + j X<-^ 


(29) 


Similarly in Region III a: > ^ 


= V., = t(fc) ( :r>e 


(30) 


For fc ^ 0 we can define two independent solutions of Eqs. (3) by, for example. 




E + m 
—ik 


Akx 


X —oo 


(31) 
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oo 
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which represent purely incoming particles from the left and right respectively. By taking appropriate linear combina¬ 
tions of and normalising we can choose two new independent solutions of Eqs. (3) 


with the properties 




(33) 


9i{-0=0 ff2(-0 = l /i(-e) = l f2{-0=0 (34) 

Note that the solutions ipi and -02 which satisfy Eq. (jsj) are everywhere real provided that k is real. We can then 
express our solution 0® in terms of a linear combination of tpi and 02 for all x and in particular in Region II \x\ < 0 


0 " = 6 


h{x) 

9i(.x) 





(35) 


We can evaluate the Wronskian of the solutions 0i,02 is constant at the point x = —^ to give 


W"[01,02] = w{-0 = /i(-e)g2(-0) - h{- 09{-0 = 1 

thus confirming that the solutions 0 i ,02 are independent for k 0. 

The wave function 0®(x) must be continuous at a; = —0 and x = ^. The overlap between Regions I and II and 


between II and III then give the following boundary conditions: 

(E-hm)(e-*'=«-hZ(fc)e*'=«) = b (36a) 

- Z(fc)e*'=«) = c (36b) 

(E + m)tik)E^^ = bhiO + c/ 2 (e) (36c) 

-ikt{k)E^^ = 6gi(0)+cg2(O (36d) 

For simplicity, write ai = /i(^) and Pi = 5i(0), so that the last two equations become 

{E + m)t{k)e‘^^^ = bai + ca 2 (37a) 

-ikt{k)E’^^ = bPi+cP2 (37b) 


Note that b and c are dependent on k as are ai and Pi. Eliminating t, b and c then re-arranging to solve for I gives 


m = 


Similarly t can be found to be 


/ fc^a 2 + {E + m^Pi + ik{E + m){ai — P 2 ) \ - 2 ik^ 
V kPa 2 — [E + my Pi — ik{E + m){ai + P 2 )) 


t{k) = 


—2ik{E + m){aiP2 — a2Pi) 
kPa 2 — {E + my Pi — ik{E + m){ai + P 2 ) 


e-2rK 


We can then use the relation 


(38) 


(39) 


W{C) = fi{C)92{i) - f2{09i{0 = “1/32 - OL 2 P 1 = 1 (40) 

to^simplify Eq (|^. It is a straightforward exercise to verify that Eqs. (|3|, ^) satisfy the unitarity condition Eq. 

We can now discuss the limit as A: ^ 0. It is apparent from Eqs. (^, ^) that provided /3i(0) y 0 the limit E ^ m, 
k ^ 0 gives 

;(0)=r(0) = -l t(0) = 0 (41) 

so that the reflection coefficients T(0) = i?(0) = 1 and the transmission coefficient r(0) = 0 . These results in the 
general case agree with those for the Schrodinger equation |^, Q. 

Using Eqs. (33a,33b) it can be seen that 5(0) = c(0) = 0 in the k = 0 limit and therefore from Eq.(|3^) the wave 
function vanishes identically for all x. Thus the only physical solution of the Dirac equation (1) for fc = 0 is the 
solution 


0(a;,A: = O) = O (42) 

unless the potential has special properties which we investigate in the next section. 

It should also be noted that as fi{x) and 9 i{x) are real at a: = ^, the quantities a^fc) and Pi{k) are also real. Hence 
from Eq. (^^ as A: —> 0, t{k) is pure imaginary as it approaches zero in agreement with the Levinson theorem Q for 
the Dirac equation provided /3i(0) y 0. 
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B. The special case /3i(0) = 0 


If we return to Equations (33b,34b) we see that as fc ^ 0 we must have 

c(0) = 0 6(0)/3i(0) + c(0)/32(0) = 0 


so 


6(0)/3i(0) = 0 

Furthermore since c(0) = 0, we must have from Eq. (34a) 

2 mt(0) = 5(0)ai(0) 


( 43 ) 


(44) 


When b{0) = 0 as well as c(0) = 0 we obtain the general case already discussed. If /3i(0) = 0, however, then as we 
approach the limit E ^ m, k ^ 0 the reflection amplitude l{k) does not satisfy 1(0) = —1 and so the wave function 
fc = 0) ^ 0. In this case we therefore have non-trivial solutions of the Dirac equation at fc = 0. This implies that 
transmission coefficient t(0) will be non-zero in this limit as will ai(0). 

For k ^ 0 with /3i(0) = 0 we can write f3i{k) = kf3i{0) . So from Eq. 


1 (0) = lim 


/32 — ai + 2ml3i{0)i 


k^o (32 + 0 . 1 + 2m/3)(0)i 


(45) 


As k is arbitrarily small (and not actually equal to zero), the Wronskian, W = ai(32 = 1-1- 0{k) so (32 = l/<ai + 0{k) 
and in the limit k 0 we have 


1 ( 0 ) = 


1 — ai(0) -I- 2miai{0)(3i{0) 
1 -I- of (0) — 2miai{0)(3'i{0) 


(46) 


We know however from Eq. (W) that 1(0) must be real. From Eq.(|46|) this means that either /3i(0) = 0 or ai(0) = 0. 
But since we are considering the non-trivial case where t(>{x, k = 0) ^ 0 (and hence we expect that t(0) ^ 0) we do 
not want ai(0) = 0 since from Eq. (H) that would imply that t(0) = 0. Thus we would like to be able to show that 


/3i(0)=0 


(47) 


and (3i(k) = 0{k^). This is not difficult to demonstrate using an argument of Lin’s [^: the wave function i/)iof 
Eq. (pd]) is a solution of Eqs. (3) subject to the fc-independent boundary conditions given by Eq. ^). So its lower 
component 


gi{x,k) = gi{x,E) 

since the Dirac equation (3) involves E explicitly not k. It follows that 

l3i=PiiE) 

which requires (3i to be an even function of k and in particular that a.s E = \Jrn3 -|- 

d(3i{k) d(3i{E) dE k d(3i{E) 

dk dE dk E dE 

at fc = 0 in agreement with Eq. (^^. 

This gives the final result for the reflection amplitudes in the special case when /3i(0) = 0: 

l-af(O) 


(48) 


(49) 


(50) 


l(0 )=-r(0) = 


1 -I- Q!i(0) 


and for the corresponding transmission amplitude from Eq. (^): 

, 2°i(0) 

l + o?(0) 

These results agree with those obtained by Clemence [||. 


(51) 


(52) 
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C. Half-Bound State 


We will now show that if the potential were to support a bound state in the limit E = m then /3i(0) = 0 so the 
scattering wave function will not vanish in the limit k ^ 0. For an asymmetric potential the following bound state 
wave function is appropriate for |a;| > 


Region I 


= s 


E + m 

— K 




X < 


Region III 


V''' = s' 




x>i 


(53) 


If the potential is such that the wave function possesses a well-defined non-zero limit as i? ^ m, k —> 0, then the 
wave function for \x\ > ^ in this limit is just proportional to 



(54) 


albeit with different constants of proportionality s, s' on the left and right. It is clear that a wave function of this form 
is non-normalisable and forms part of the continuum. 

The scattering solutions which tend to the solutions Eq. (Q) in the zero-momentum limit will therefore have a 
lower component which vanishes for sufficiently large |a;|. From Eq. ( ^^ this implies that at x = ^ 

6(0)/3i(0)+c(0)/32(0) = 0 (55) 

while at X = —^ using Eq. (^) we have c(0) = 0. Since ip^{k = 0) is not zero for a half-bound state, 6(0) ^ 0 and 
hence 


/3i(0)=0 


(56) 


An example of a half bound state in an asymmetric potential is given in the Appendix together with an explicit 
demonstration that Pi{k) is of order for small k when the condition /3i(0) = 0 holds . 


D. Symmetric Potentials 

When the potential is symmetric so that V{x) = V{—x) we can find more stringent conditions on 1(0), r(0) and 
t(0). In the two-component approach, the behaviour of the wave function under the parity transformation x —> —x is 
given by: 


ip'i-x^t) = a^ip{x,t) (57) 

It follows that we can define an even wave function ip+ix) under parity as one with an even top component and an odd 
bottom component whereas an odd wave function ip-lx) has an odd top component and an even bottom component. 
The wave function ip^ for the bound state given in Eq. (|53| ) must now be either an even solution ip+ or an odd solution 
Ip-. First let us assume that it is even. 

Then in the limit of a half-bound state at E = m, (k —> 0) the solution remains even. As k ^ 0 the scattering 
solution ^p'^ will also be even. Thus from Eqs. (Q, ^ we have 

l + l{0) = t{0) (58) 


From the unitarity relation we also know that 


1 ( 0)2 -h t ( 0)2 = 1 = 1 ( 0)2 + (1 + ^( o ))2 


therefore 


So either 1(0) = 0 or 1(0) 


1(0)2-h 1(0) = 0 


(59) 


— 1 in agreement with Eq.(^). We know that 1(0) ^ 


— 1 as {k = 0) 7 ^ 0. Hence 
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;(o) = 0 


(60) 


and the transmission coefficient 


r(0) = 1 (61) 

Using Eq ( ^^ we see that for an even half-bound state we must have ai(0) = 1 while for an odd half bound state we 
have ai(0) = —1. 

So we obtain the result that when a symmetric potential supports a half bound state, a transmission resonance 
T = 1 occurs for an incident particle with arbitrarily small momentum. This agrees with our previous result for 
reflectionless scattering by a repulsive potential V(x) where its attractive counterpart U(x) = — V(x) is supercritical 
[Q, that is to say U{x) has a half-bound state ai E = —m. To see this note that Eqs. (3) are invariant under the 
(charge conjugation) transformation 


E^-E V^-V f^g g^f 

so it follows that V{x) has a half-bound state at E = m when U{x) has a half-bound state at E = —m. 


(62) 


DISCUSSION 

We have now generalised the results for scattering in one dimension in the Schrodinger equation to the Dirac 
equation as we intended. But we are physicists not mathematicians: consequently our results are not yet as complete 
as those proven for the Schrodinger equation. Clemence |^, however, has shown that the class of potentials for which 
our results are true in the Dirac equation can be extended to include potentials which do not vanish for \x\ > His 
results require the potentials V(x) to satisfy 


/ (1 -I- |a;|)|U(a:)|dx < oo (63) 

J —OO 

As stated in the Introduction a half-bound state at E = m can arise in two ways in the Dirac equation. These can 
most easily be distinguished by the examples of an attractive well for which V (x) < 0 and a repulsive barrier for which 
V(x) > 0, although it may be more difficult to characterise which is which for a complicated potential. In the case of 
an attractive potential a half-bound state with E = m corresponds to a non-relativistic zero energy resonance. For 
example in the case of a square well V(x) = — Vq, |a;| < a, U(a:) = 0 elsewhere one occurs at the threshold for the first 
odd state Vq = /2ma?. In the case of a repulsive potential a half-bound state occurs as we have just seen when the 

corresponding attractive potential U{x) = —V{x) is supercritical. For the square barrier V{x) = Vq, |a:| < a,V{x) = 0 
elsewhere supercriticality first occurs when Vq = m -I- -I- Tr^/da^ Note that Vq > 2m before supercriticality 

can occur. 

Over 70 years ago Klein 0 discovered that a Dirac particle could tunnel through a potential barrier V with 
V > 2m. In this paper we have confirmed that tunnelfing will always occur in the Dirac equation if a potential 
barrier V{x) of short range is strong enough so that U{x) = —V{x) is supercritical. The generic phenomenon whereby 
fermions can tunnel through barriers without exponential suppression we have called “Klein Tunnelling” [^ . Even 
strong long range repulsive potentials in the Dirac equation seem to have this property: in three dimensions Hall and 
one of us (ND) 0 have shown that Klein tunnelling is also associated with supercriticality for Coulomb potentials. 

Acknowledgments We would like to thank Peter Bushell and Alex Sobolev for their help. 


I. APPENDIX 

In order to illustrate scattering off an asymmetric potential we shall consider one of the few examples which can be 
solved analytically. We shall use a double delta potential barrier which comprises two unequal Dirac delta functions: 

V{x) = Xd{x) + fi d{x — a) (64) 


where \ ^ g. and X, g > 0. 
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A. Scattering Coefficients 


The wave function for a: < 0 is 


while for 0 < X < a it is 


and for x > a 


'ipix) = 


E + m 


E + m 
—ik 


ikx 

E + m\ 

ik J 

,^kx p 

/ E + m 
ik 

E + m 

1 a^kx 

—ik 

r 


The discontinuity condition on ip{x) at the first barrier at x = 0± is 

The second discontinuity condition at x = a± is derived by replacing 0± with a± and A with /i. 

The reflection and transmission amplitudes I and t can then be calculated to give: 

imk{cos /i sin A + e^“^cos A sin fi) + — l)sin A sin fi 


and 


I = -- 


t = 


^ 2 ^g 2 iafe _ i^gin A sin ^ + k?cos[X + /r) + iEk sin{\ + 


j^ 2 ^g 2 iafe _ \ giji ^ _|_ k^cos{X + ^) + iEk sin{X + /i) 

Using E = we can write for small k 

—im{sin{X + /i) + 2ma sin X sin fi) + 2anik{am sin X sin /r + cos X sin /r) + 0{k^) 


I = 


t = 


im{sin(X + /x) + 2ma sin X sin fi) + k(cos{X + fj.) — 2a‘^m?sin X sin fj) + 0{k^) 

k 


im(sin{X + fi) + 2ma sin X sin fi) + k{cos{X + fi) — 2a?nT?sin X sin fi) + O^k?) 
From Eqs. (^, |^) it is easy to see that in general as fc —> 0 


in agreement with Eq. (^). If however 

sin{X + fi) + 2ma sin X sin fi= 0 


then 


and 


I 


am sin{X — fi) 
cos(X + fi) + am sin{X + fi) 


1 

cos(A + fi) — 2a?m‘^ sin X sin fi 


It is easy to show that I and t given above do indeed satisfy 

If+ |t|^ = l 


provided that sin{X + /i) + 2ma sin X sin fi= 0. 









B. Exceptional Case 


The exceptional case in the proof above occurs when /3i(0) = 0. We shall therefore calculate ai(0) and /3i(0) for the 
double delta potential. From Eq. (|^) we consider the solution of the Dirac equation which takes the values ^ q ^ at 
X = The wave function for a; < 0 thus has the form 


{E + m)ip{x) 


(E + to) cos k(x + 
k sin k(x + 


(77) 


while for 0 < X < a it is 


. ^ \ \ f (E + to) cos kx\ , e / (E + m) sin kx \ 

(E + to)^(x)= 7(^' fcsanfcx j+% -kcoskx ) 


and for x > a we can write 


/ rn , \ I / \ f (E + m) cos k(x — 0 \ . f (E + m) sin k(x — $) 

' /V / 1 ksmkyx — ^) j \ —kcoskyx — ^) 

So at X = ^ we see that 

ai{k) = a Pi{k) =—kr/{E + m) 

For small k we calculate from the discontinuity conditions that 

cr = [cos(A + /i) + 2 to^ sin{X + /a) — 2ma sin fi cos X + 4aTO^(^ — a)sin ^ sin A] + 0{k^) 


(78) 


(79) 


(80) 

(81) 


and 


kr = 2m [(sm(A + /r) + 2am sin X sin p,)] + 0{k^) (82) 

Note that neither a nor kr has any term of order k. As A: —> 0 we obtain 

/3i (0) = — [(szn(A + fi) + 2am sin X sin p,)] (0) = 0 

so the exceptional case given by Eq. (^ above indeed satishes /3i(0) = 0. Furthermore when /3i(0) = 0 it is easy to 
see that 


ai (0) = cos(A + p.) — 2am sin ^ (cos A + 2am sin A) 
From Eq. (j^) above the transmission coefficient in the exceptional case when 

/?! (0) = — [(sw(A + ^) + 2am sin X sin fi)] = 0 


can be expressed in terms of ai (0) 


2ai(0) 

1 + ai(0) 


After some tedious manipulation we find that 


1 + [ai(0)]^ = 2(1 + 2amsinXcos X + 2a^m^ sin^X) 

= 2(cos(A + /i) + 2am sin X cos /i) (cos(A + /i) — 2a?'m? sin X sin /r) 
= 2ai(0)(cos(A + /i) — 2a'^m? sin X sin fi) 


(83) 


(84) 


(85) 


So 


t = 


2ai(0) 


2ai (0) [cos(A + ^) — 2a'^m‘^ sin X sin 
in agreement with Eq. 76). Similarly it can be shown that Eq. ([^ for the reflection coefficient agrees with Eq. 


( 86 ) 
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C. Bound States 


Let us now consider the asymmetric potential well 

U(x) = —V(x) = —\6{x) — fi S{x — a) 
This will have bound states with a wave function for a; < 0 of the form 


while for 0 < x < a 


and for x > a 


ip{x) = 7 


ip(x) = 


— K 

m — E 


m — E 


m — E 


ip(x) = s 

The discontinuity condition for the first delta well is 


m — E 




(87) 

( 88 ) 

(89) 

(90) 

(91) 


Note that this differs from the condition for barriers only in that A ^ —A. The second discontinuity condition follows 
with 0± ^ a± and A —> /r. This leads to the following four equations: 


— i ?)(7 -I- 6)sin X = —K 

(92a) 

- «:(—7 -1- 6)sin X = ni — E 

(92b) 

+ Se-^'^)sin fi = SKc-^'^ 

(92c) 

,a« ^e-“«)szn fi = s{m- E)e-^^ 

(92d) 


7 and 5 can be found from the first two equations (88a, 886) to be: 

Esin A — Kcos A 


7 = 


S = -- 


:A 


(93) 

A. 

Eliminating s from (88c, 88d) leads to: 

(kcos ^ — Esin ^) + Sm sin ^ = 0 (94) 

We thus obtain: 

e^°‘^{Kcos A — Esin X){kcos /i — Esin /r) — m^sin A sin /r = 0 (95) 

Re-arranging gives 

, K^e^°^'^cosXcos u + {e^°‘'^E'^ — ni?)sin\sin u 
sin{X + /i) =-- (96) 

At supercriticality E = —ni, k = 0 giving 

sin{X + ^) + 2ma sin X sin fx = 0 (97) 

in agreement with the exceptional condition /3i(0) = 0. 

When X = fj, we obtain a symmetric potential. If sin{X + fi) + 2ma sin X sin /i = 0 then either sin A = 0 and 
q:i( 0) = 1 or tan A = —1/ma and ai(0) = —1. In both cases the transmission coefficient T = 1 in agreement with our 
previous result Q for supercritical symmetric potentials. 


11 






[1] Fadeev, L.D.,Trudy Mat. Inst. Stekl. 73 314(1964) [Amer. Math. Soc. Translation 2 139 (1964)] 

[2] Bohm, D., Quantum Mechanics, (Prentice Hall, Englewood Cliffs, NJ, 1951), p. 286. 

[3] Bolle D., Gesztesy F. and Wilk S.F.J., J. Operator Theory 13 3 (1985) 

[4] Sassoli de Bianchi, M., J. Math. Phys. 35 2719 (1994) 

[5] Clemence D.,Inverse Problems 5 269 (1989) 

[6] Newton, R. G., Scattering Theory of Waves and Particles, (Springer-Verlag, Berlin 1982), p.280. 

[7] Dombey, N., Kennedy, P. and Galogeracos, A., Phys. Rev. Lett. 85, 1787 (2000) 

[8] Galogeracos, A., Dombey, N. and Imagawa, K., Yadernaya Fiz. 159, 1331 (1996); Phys. At. Nuc. 159, 1275 (1996) 

[9] Gohen-Tannoudji, C., Diu, B. and Laloe, F., Quantum Mechanics, Volume 1, (New York, Wiley, 1982), p. 359-66. 

[10] Barthelemy, M.-C, Ann. Inst. Henri Poincare 7 115 (1967) 

[11] Galogero, F., Variable Phase Approach to Potential Scattering, (Academic Press, New York, 1967), p. 113. 

[12] Senn, P., Am. J. Phys. 56, 916 (1988) 

[13] Lin Q-G, Fur. Phys. J. D7, 615 (1999) 

[14] Klein O., Z.Phys. 53 157(1929) 

[15] Galogeracos A. and Dombey N., Int J Mod Phys A14 631(1999); Phys. Rep. 315 41(1999) 

[16] N Dombey and R L Hall, Phys. Lett. 474 1 (2000) 

[17] Kennedy,P., J.Phys A35, 689 (2002) 


12 



